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Abstract 

Using the modified factorization method employed in ||, we construct a new 
class of radial potentials whose spectrum for I = coincides exactly with that of 
the hydrogen atom. A limiting case of our family coincides with the potentials 
previously derived by Abraham and Moses |6|. 
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Ph! 1 Introduction. 



In almost all exactly soluble spectral problems of the energy operator, the possibility 



of obtaining an exact expression for the energy levels and the corresponding eigenstates 
stems from the group theoretical symmetries which lead to the algebraic factorization 
method This method, until quite recently, appeared to have been fully explored. 

^ ■ Yet, from time to time, some new potential arises for which the factorization method can 

be applied in a new way ||]. Recently, a variant of the factorization method proposed 
by Mielnik || has permitted a description of a new class of potentials whose spectra are 
identical to those of the harmonic oscillator Below we apply this method to show 

that there is a one-parameter class of radial potentials, different from a Coulomb field for 
which the sequence of energy eigenvalues for / = is exactly as for the hydrogen atom. 

2 Hydrogen potential: standard and generalized fac- 
torizations. 

As is well known, the hydrogen atom eigenproblem: 

-^V 2 -^j^(r) = ^(r), (2.1) 
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which, after the separation of the angular variables ip(r) = Y(9, ip)R(r) leads to a sequence 
of radial problems: 



d 2 | 1(1 + 1) 



dr 



rR = XR, 



(2.2) 



where I = 0, 1, 2, . . . labels the angular momentum eigenvalues, r is a new dimensionless 
radial coordinate, and the functions R = {R(r)} (0 < r < +oo) form a Hilbert space 
H r with the scalar product defined by: (R, R') = 4tt J +o ° R(r)R'(r)r 2 dr. It is also well 
known that the radial Hamiltonians 



H, = - 



^<f_ 1(1 + 1) 2 

dr 2 r 2 r 



(2.3) 



admit a sequence of 'factorized forms' which allows us to find very simply the energy 
levels iTH. These forms are 



Hi = afai — — 
Hi = ai + ia l+1 - 



(i + iy 



where 



ai 



d I 1 

dr r I 



d I 1 
dr r I 



r. 



(2.4) 
(2.5) 

(2.6) 
(2.7) 



Now, following the method in H], we would like to ask two questions. Is this representation 
unique? Are ai and a~l the only first-order differential operators for which ( |2.4]) - (|2.5| ) hold? 
It turns out that the answer is negative. Put: 



A, 



At 



dr 



and demand that the formula analogous to 



AfAi = H, 



should be again valid: 



1 



This leads to 



#(r)+tf(r) 



I 2 ' 
1(1 + 1) 



2 1 
- + 72- 



Taking into account that we already have one particular solution, (3i = (l/r) 
general solution of (|2.11|) can be easily obtained. Denote 



01 



I 1 

r I 

2 



<Pi{r)- 



(2.8) 
(2.9) 

(2.10) 

(2.11) 
(1//) the 

(2.12) 



Then 



21 2 

if 1 J i 2 

-<Pl + —<Pi ~ j<Pi + <?i 
r I 



0. 



After introducing a new function, Xi = 1/0/, we obtain: 



, 21 2 
X + -X t - -Xi + 1 = 
r / 



with the general solution 
and, therefore, 



ll 



y 2l e- 2 yl l dy)r- 2l e 2r l l - 



AW 



/ 1 

- T + 



r 2l e -2r/l 



r I li — Jo V 2le 2y l l dy 
The commutator of the new operators Ai and Af is not a number: 

[A h Al]=2Pl{r). 

Thus, we can apply the factorization method in a new way. 

Taking the product Ai Af , instead of Af Ai in ( [2.10|) , one has 

AAf = Af Ai + [A h Af]=H l + ± + 2#(r) = Hi-i + ± 

where Hi-i is a new Hamiltonian: 



with 



If 



VU(r) 



= H l + 2#(r) 



2 1(1-1) d 

+ :T 

r r ar 



dr* 



2 r 2l e -2r/l 



H ~ JoV 2l e 2y/l dy 



, 1 = 1,2,3, 



2Z+1 



7< > (2/)! 



(2.13) 



(2.14) 



(2.15) 



(2.16) 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



or 7; < for a fixed I, the third term has no singularities. Furthermore, V5_i(r) — > as 
r — > oo, and so we obtain a one-parameter family of new self-adjoint Hamiltonians. To 
find the spectra of Hi-i note that: 



H^Ai = ( A t Af - — ) Ai = Ai[ AfAi - — 



(2.21) 



This implies that if {R n i} are eigenvectors of Hi with eigenvalues A n , {AiR n i} are eigen- 
functions of Hi_\ with the same eigenvalues. Furthermore 



AlRnl = aiRni + 



r 2l e -2r/l 



li-Ioy 2l e- 2y/l dy 



R 



Til 



(2.22) 



are square-integrable functions. They are orthogonal due to 

(AiRnijAiRn'i) = (AfAiRni, Rn'l) = {^n + j^j <W ■ (2.23) 

As one can also check, the operator A\ maps the continuous spectrum subspace of 
Hi into the corresponding continuous spectrum space of However, the operator A\ 

does not map the Hilbert space 7i r of radial wave functions into the whole of 7i r . What 
remains to be examined, similarly like in ||, is the 'missing vector' Ri /_i orthogonal to 
all vectors of form AiR(R e Ti r ): 



(Ri i-uAR) = (A+Ri i-i,R)=0. 

R 



(2.24) 



This means that Ri ;_i is obtained form the first-order differential equation 



AfRi i-i 



1 



-£+ AW 



rRi i-i = 



whose solution is 



Ri i-i 



C\r l 1 e r l l 



(2.25) 



(2.26) 



H ~ Ioy 2l e 2y / l dy' 

One immediately checks that Ri ;_i is an eigenvector of Hi_i with the eigenvalue —(l/l 2 



1 



Hi-iRi i-x = AiA, — — ) Ri i-i = — -rR 



For any 



21+1 



1 ~ 



I 2 ' 



H l-l- 



(2.27) 



li > (21)1 



or 7; < 0, 



equation (|2.26|) is a square-integrable function defining the new ground state of the mod- 
ified potential. Hence, one can see that when 



7/ > (2/)! (~ 



21+1 



or 7; < 0, 



Hi-i is a new family of self-adjoint radial Hamiltonians with the same discrete spectrum 
as Hi-\. 

The limiting case 7, = (2Z)!(//2) 2m is worth attention. When 7, = (2Z)!(//2) 2m the 
third term in (|2.16|) tends to a constant when r tends to 00. Hence, AiR n i is remaining 
square-integrable and VJ_i(r) tends to zero when r tends to 00. However, the function 
defining the new ground state is no longer square-integrable. Therefore, in this case we 
obtain a new potential in which the lowest energy level is missing comparing with Hi_\. 



4 



3 The case / = 1. 



By taking I = 1 (i.e., starting our procedure from the conventional Coulomb potential with 
the simplest centrifugal term) we pass to new radial functions different from the Coulomb 
potential but for which the sequence of energy levels corresponding to the vanishing 
angular momentum is exactly the same as for the Hydrogen atom. Our new potentials 
still depend of one arbitrary parameter: 

71 - 1/4 + (r 2 /2 + r/2 + l/4) e - 2r J ' ^' ' 

It is interesting to notice that for 71 > 1/4 the potentials ( |3.1|) have the same singularity 
at r = and the same asymptotic behaviour at r — > 00 as the Coulomb field (see Figure 
1). The corresponding ground states are plotted on Figure 2. 
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Figure 1: The potentials Vo(r) for various values of 71. 



The critical value of 71 — > 1/4 is of interest. Our potential (|3.1| ) then tends to: 



16r(r + 1) 
(2r 2 + 2r + 1 



|2' 



(3.2) 



All the eigenvectors for 71 — ► 1/4 converge to some square- integrable function, except the 
first one ( |2.26| ) (missing vector) which then becomes non-normalizable. Therefore, the 
potential (|3.2| ) is missing the lowest energy level comparing to the hydrogen atom. As 
immediately seen, this case reproduces the one previously found by Abraham and Moses 
by applying the Gelfand-Levitan method (see 0). For 7 7^ 1/4, the radial potentials 
(|3.1|), as far as we know, have not been known before. 
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Figure 2: The ground states i?io(r) for various values of 71. 
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